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When a colloidal particle adheres to a fluid membrane, it induces elastic deformations in the 
membrane which oppose its own binding. The structural and energetic aspects of this balance are 
theoretically studied within the framework of a Helfrich Hamiltonian. Based on the full nonlinear 
shape equations for the membrane profile, a line of continuous binding transitions and a second line 
of discontinuous envelopment transitions are found, which meet at an unusual triple point. The 
regime of low tension is studied analytically using a small gradient expansion, while in the limit 
of large tension scaling arguments are derived which quantify the asymptotic behavior of phase 
boundary, degree of wrapping, and energy barrier. The maturation of animal viruses by budding is 
discussed as a biological example of such colloid-membrane interaction events. 
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I. INTRODUCTION 

Arguably the most important structural component of 
all living cells is the phospholipid bilayer. It combines 
the two diametrical tasks of partitioning - thereby orga- 
nizing the complex hierarchy of intracellular biochemical 
environments - while at the same time providing con- 
trolled transport mechanisms between neighboring com- 
partments The size of particles being transported 
spans several orders of magnitude, ranging all the way 
from sub-nanometer ions, whose passage through the bi- 
layer is facilitated by protein channels, up to micron-sized 
objects engulfed by the large scale membrane deforma- 
tions occurring during phagocytosis. 

Since cell survival depends on a meticulous balance of 
these processes, they are actively controlled and require 
metabolic energy to proceed. Still, there are cases where 
they happen passively as a result of generic physical in- 
teractions, for instance, a sufficiently strong adhesion be- 
tween the particle about to be transported and the mem- 
brane. An example that has been studied in extensive de- 
tail is the route along which many animal viruses leave 
their host cell Q- After entering the cell (typically via 
receptor mediated endocytosis or other active processes 
0) and completion of the viral replication steps, the new 
virions have to get out again. Many virus families accom- 
plish this by the plasma membrane wrapping their nu- 
cleoprotein capsid and pinching off ( "budding" ) — a step 
which not only sets them free but by which they also 
acquire their final coating. A particularly clean model 
case is provided by Semliki Forest Virus, in which case 
the binding between capsid and membrane is promoted 
by viral ("spike"-) proteins 0,|E@l- 

A different realization of such a wrapping event is pre- 
sumably exploited by an efficient gene transfection sys- 
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tern proposed a few years ago by Boussif et al. Q- There, 
DNA is complexed by the cationic branched polymer 
polycthylcniminc into a globular complex, which then en- 
ters the cell (as deduced from reporter gene expression). 
Rather than particular targeting signals, the slight net 
positive charge of the complex is believed to trigger ad- 
hesion and membrane penetration via an electrostatic in- 
teraction with negatively charged regions of the plasma 
membrane. 

Finally, a great deal of biophysical experimental tech- 
niques involve the attachment of microbeads to mem- 
branes. For instance, one way to measure cellular 
tensions involves pulling a thin tether with an optical 
tweezer which grabs a bead adhering to the membrane 
0. Cell membranes and subcellular organelles are rou- 
tinely probed with an atomic force microscope 9]. And 
a classical experiment on surface dynamics and locomo- 
tion of cells involves monitoring the centripetal motion 
of a surface adherent particle In all cases the object 
adhering to the cell will locally deform its plasma mem- 
brane, which can be crucial for interpreting the experi- 
mental results. For example, in the case of centripetal 
bead motion it has been noted that adhering particles 
may actually become engulfed by the cell without involve- 
ment of the endocytosis machinery, if only the membrane 
tension is low enough [lT| . 

Wrapping and budding processes occur very frequently 
in cells, but even though they are extensively studied ex- 
perimentally, the high complexity of the real biological 
situation renders a clear extraction of underlying phys- 
ical principles very difficult; it is not even obvious or 
undisputed whether an explanation in terms of such prin- 
ciples is possible. A better look at the physics is there- 
fore provided by more easily controllable experiments on 
the adsorption of colloids onto model lipid bilayers. The 
degree of wrapping of a colloid by a giant phospholipid 
vesicle has for instance been studied by Dietrich et al. 

. These authors showed that it can be quantitatively 
understood in terms of a balance between adhesion and 
elastic energy. Colloids will change the shape of the vesi- 
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cle they adhere to, and this may give rise to attractive 
interactions between them |l3j . If the membrane tension 
is low enough or the colloid sufficiently small, bending 
of the membrane will become an important contribution 
which can be strong enough to completely suppress adhe- 
sion 0, 0, 0] . The bending stiffness will also prevent 
any "kink" in the membrane profile at the line of con- 
tact, such that the notion of a contact angle only remains 
meaningful in an asymptotic sense and is replaced by the 
concept of contact curvature 0, 0, 0, |2(j . For vesicles 
adhering on flat surfaces, this can be measured by visu- 
alizing the contact zone via reflection interference con- 
trast microscopy , but for adhesion of membranes on 
strongly curved substrates - in particular, small colloids 
- this option is not available. Here, computer simulations 
can provide means for a close-up study [22j, but to date 
it is still difficult to perform a quantitative analysis relat- 
ing the observed geometry to the underlying mesoscopic 
elastic properties. 

Instead of directly measuring the local deformations, 
one can study their indirect consequences and from there 
attempt to deduce some of their properties. This of 
course requires theoretical modeling for bridging the gap. 
For instance, lateral membrane tensions can be inferred 
from the force required to pull a thin tether from the 
membrane . Recent theoretical work by Derenyi et al. 
shows that the force as a function of tether length is quite 
subtle and shows an initial oscillatory structure before 
exponentially settling down to the asymptotic value |23| . 
Very useful information might be obtained from this, for 
example an independent estimate of the tether width. For 
less strongly bound beads Boulbitch has related the un- 
binding force to the adhesion energy, membrane bending 
stiffness and the elasticity of an underlying elastic net- 
work (like the cytoskeleton) by describing the associated 
membrane deformation within a small gradient approx- 
imation [24|. And for the case of colloids adhering to 
quasispherical vesicles a simple ansatz for the membrane 
shape has been shown to yield a structural phase diagram 
indicating when the colloid is free, partially wrapped and 
fully enveloped (ltsf. 

This paper extends the work of Ref. [25( in develop- 
ing a detailed theory of the local wrapping behavior in 
the case of a constant prescribed lateral membrane ten- 
sion. This "ensemble" is particularly relevant for the 
biological situation, since almost all cells constantly ad- 
just the amount of lipids in their plasma membrane in 
order to maintain its lateral tension at some specific set 
point [2(| . The situation to be discussed here is similar 
to the one treated by Boulbitch [2_^|, but it will neither 
be restricted to small membrane deformations nor to a 
two-dimensional modeling of the geometry (which would 
be appropriate for long cylindrical colloids). After set- 
ting up the problem and identifying the relevant ener- 
gies in Sec. [H] the full nonlinear shape equations of the 
membrane profile are studied in Sec. IIIII leading to the 
structural wrapping phase diagram. Sec. HVI then treats 
the small gradient expansion (i. e., the linearized theory), 




FIG. 1: Illustration of the wrapping geometry and membrane 
parameterization. A membrane adheres partially to a spher- 
ical colloid of radius a with a degree of wrapping given by 
z = 1 — cos a. Cylindrical symmetry around the eh axis is as- 
sumed throughout the paper. The possibility of "overhangs" 
requires going beyond a Monge-like parameterization, which 
would gives the height h as a function of radial distance r. 
The choice followed in this work is to specify the angle ijj 
with respect to the horizontal as a function of arc-length s. 

identifies its range of validity as the regime of low mem- 
brane tension, and derives various asymptotically exact 
results. For the opposite regime of large tension several 
scaling predictions are deduced in Sec. |V] and validated 
against the numerical results from Sec. Ifffl In the final 
Section lvTl the results obtained throughout the paper are 
discussed in the framework of the biological example of 
virus budding. 

II. GENERAL ENERGY CONSIDERATIONS 

The aim of this paper is to describe the local deforma- 
tions of a flat fluid membrane upon adsorption of a spher- 
ical colloid, as depicted schematically in Fig.^ It will be 
assumed that this process can be understood as a balance 
of the following three energy contributions: (i) Adhesion 
is driven by a contact energy per unit area, w, and op- 
posed by (ii) the requirement to bend the membrane as 
well as {Hi) the work of pulling excess membrane toward 
the wrapping site against a prescribed lateral tension, a. 
Since the description will not aim at a microscopic under- 
standing, continuum elasticity theory is taken as a basis; 
in particular, the bending energy per unit area will be 
described using the standard Helfrich expression [27j 

ebend = ^k(ci + c 2 - c ) 2 + KCiC 2 , (1) 

where c\ and c 2 are the local principal curvatures of the 
two-dimensional membrane surface |28| , Co is the sponta- 
neous curvature of the membrane and k and R are elastic 
constants (with units of energy). In the following a sym- 
metric membrane is assumed, i. e. Cq = 0; and since no 
topological changes will be considered, the second term in 
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Eqn. (JIJ can be dropped [2^| . The tension energy is given 
by the lateral tension, a, times the excess area. Note that 
from tension and bending constant one can construct a 
length A - specific to the membrane - according to 

Membrane deformations on a length scale smaller than A 
predominantly cost bending energy, while deformations 
on a larger scale pay mostly in tension. An interesting 
situation occurs when the two intrinsic lengths of the 
problem - the colloid radius a and the membrane length 
A - are of the same order, because then bending and ten- 
sion contributions are comparable. In fact, for a typical 
cellular membrane tension of a ~ 0.02dyn/cm |2g and 
a typical bending modulus of k ~ 20 k^T (where k^T is 
the thermal energy) one obtains A ~ 64 nm. Viral capsids 
are about this big |2(, therefore this biological situation 
sits squarely in the crossover regime. 

Assuming that the adhering membrane remains in its 
fluid state, its energy can be calculated easily once the 
degree of wrapping, z = 1 — cos a, is specified, because 
the shape of the membrane is known. The area of the 
colloid covered by membrane is given by A a d = 2ira 2 z, 
which gives an adhesion energy of E a< \ — — wA a d = 
—2na 2 zw. Using Eqn. l|T]l. the bending energy is found 
to be -Ebcnd = |/t(l/a + l/a) 2 A a d = Attzk. Finally, 
the work done against a lateral tension a is proportional 
through a to the excess area pulled toward the wrapping 
site, which is AA a( i = 7ra 2 z 2 , giving the tension energy 
Eton = ira 2 z 2 a. 

It turns out to be advantageous to measure energies in 
units of the bending constant k and lengths in units of 
the colloid radius a. This suggests the definition of the 
following three dimensionless variables: 

~ E _ 2wa 2 , _ era 2 
E := — , w := . and cr := , (3) 

TTK K K 

where numerical factors of tt and 2 have been introduced 
for later convenience. Note that the crossover relation 
A ~ a corresponds to a ~ 1. In terms of these re- 
duced variables the total energy of the colloid-membrane- 
complex is given by 

E = -(w - 4)z + az 2 + E {lee , (4) 

where Ef TCC — E{ vcc /ttk is the (dimensionless) energy of 
the free part of the membrane. This term is not so easy 
to calculate, since the membrane shape is not known a 
priori. Even though the main purpose of this paper is to 
determine £f rGe and understand its implications on the 
wrapping process, it still proves instructive to have a brief 
look at the problem while ignoring Ef ree . Minimizing 
Eqn. Q with respect to z (and noting that the minimum 
may also be located at the boundaries z — or z = 2) 
it is found that colloids do not adhere if w < 4, since 
they cannot pay the bending price. Once w > 4, they 
start to adhere by first being partially wrapped. Full 



envelopment occurs only if w > 4 + Ad. In between, the 
degree of partial wrapping is z — (w — 4) /2a, and both 
transitions from free to partially wrapped as well as from 
partially wrapped to fully enveloped are continuous. 

There is one case in which the simplification Ef rcc = 
in fact holds rig orously, and that is the case of zero ten- 
sion, a — |14L Il5j. This can be seen as follows: Con- 
sider the two-parameter family of cylindrically symmetric 
surfaces r(h) = a\ cosh[(h/a — a2)/flsi], which are called 
catenoids. It is easily verified that the particular choice 
a± = z(2 — z) and 02 = z — 1 + z{ ~ 2 2 z ~ > In ^z-j implies that 
this surface smoothly touches the colloid with a degree 
of wrapping equal to z. The important point is that such 
catenoids are minimal surfaces which have zero mean 
curvature at every point |28| . and hence no bending en- 
ergy. Since (i) for a — the only possible energy of the 
free membrane comes from bending, and (it) the bend- 
ing energy is positive definite, the catenoid is in fact the 
minimum energy shape. The above simplified descrip- 
tion of the problem is thus exact in the case a = 0, i.e., 
for w < 4 the colloid is not wrapped, and for w > 4 it 
is completely wrap ped, with no energy barrier impeding 
the transition |30| . 



III. THE FULL SOLUTION OF THE SHAPE 
PROFILE 

By setting Ef ICC = one neglects in particular any po- 
tential effect that the strongly curved part of the mem- 
brane close to the line of contact may have on the wrap- 
ping behavior. As a next step it is therefore tempting to 
approximate -Ef ree by a phenomenological line energy |30| . 
However, neither the relation between the line tension 
constant and the membrane properties k and a would be 
known, nor is the implied dependency on the degree of 
wrapping - namely ~ sin a — y/ z(2 — z) - supported by 
more careful studies. (On a scaling level and in the high 
tension limit this question is revisited in Sec. WV) 

In order to draw sound conclusions, an explicit treat- 
ment of the membrane deformation is needed. Three 
approaches are discussed in the following: In the present 
section the equilibrium membrane shape is determined 
by minimizing its energy. Since this leads to complicated 
nonlinear differential equations, one largely has to rely 
on numerical solutions here. In the second approach, 
discussed in Sec. IIVI a restriction to small membrane 
deformations renders these differential equations linear, 
and they can be solved exactly then. However, the lat- 
ter approach per construction is limited to fairly small 
degrees of wrapping. Hence, its range of validity is not 
obvious and has to be checked against the nonlinear re- 
sults. Finally, in Sec. 0a scaling analysis is used to shed 
light onto the high tension regime. 
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A. Energy functional and shape equations 



tions: 



The energy of the free membrane is the surface in- 
tegral over the local bending and tension contributions 
and is thus a functional of the shape. Using the cylindri- 
cally symmetric angle-arc-length parameterization from 
Fig. ^ the two principal curvatures are found to be 
(sin ip)/r and ip, where the dot indicates a derivative with 
respect to the arc- length s. The energy functional can 
then be written as 



-Efrec = / ds L(ip,tjj,r,r,h,X r ,Xh) , (5) 
Jo 

where the Lagrange function L is defined by 



L 




25-/ 

H — o(l — cosip) > 
r ) a z 

+ \ r (r — cos ip) + \h (h — sin ip) . (6) 



The expression in curly brackets contains the bending 
and tension contributions, while the two additional terms 
enforce the nonholonomic parameterization constraints 
r = cos ip and h — sin ip by means of the Lagrange para- 
meter functions X r (s) and A/j(s). Measuring all lengths 
in units of a, it can be verified that Ef rcc depends para- 
metrically on a (and through the boundary conditions 
on z) but on nothing else. This is important, because 
it ensures that w and a remain the relevant axes for a 
structural phase diagram and no new independent vari- 
able is introduced by the free part of the membrane. It 
also implies that at given a and z the membrane shape 
scales with the colloid size a. 

The Lagrangian L is independent of the arc-length 
s, therefore the corresponding Hamiltonian is conserved. 
Moreover, numerically one usually integrates systems of 
first order differential equations. Both observations sug- 
gest to switch to a Hamiltonian description. After defin- 
ing the conjugate momenta 
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a Legendre transform yields 

H = ijjpj, + rp r + hp h - L 

= r (l-cos^ 

4r r a 

+ p r cos ip + ph sin ip . 



(8) 
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' — +PrJsinV', (9d) 

^j- (1 — cos ip) , and (9e) 
(9f) 



Shape equations of this kind have been studied exten- 
sively in the past, leading, among many other things, to 
a very detailed understanding of vesicle conformations (in 
which case one also needs to fix surface and volume by 
additional Lagrange multipliers). For a detailed review 
on this subject see Rcf. 32] . 



B. Boundary conditions 

The situation to be studied is a colloid wrapped by an 
initially flat membrane. The boundary conditions thus 
have to ensure that the membrane touches the colloid 
smoothly and becomes asymptotically flat at large radial 
distances. At contact, s = 0, the following must evidently 
hold (see again Fig. [TJ: 



r(0) = a sina , 

h(0) = —a cos a , and 

tp(0) = a . 



(10a) 
(10b) 
(10c) 



The notion of asymptotic flatness can be enforced by re- 
quiring the angle tp(s) and all of its derivatives to vanish 
in the limit s — ► oo 33] . However, it suffices to demand 
this for the angle ip and the meridinal curvature tp, 



lim tp(s) = and 
lim tp(s) = , 



(11a) 
(lib) 



which implies that curvature and tension energy density 
vanish if one moves away from the site where the mem- 
brane shape is perturbed by the adhering colloid. If ip(s) 
vanishes sufficiently rapidly (as it does for a > 0, see 
Sec. HV() . all contributions beyond some large distance S 
in arc-length will be negligible. A convenient way to ex- 
ploit this is the following: Choose an upper arc-length S 
and impose the zero angle condition there. Hence, varia- 
tions of S and ip(S) are not permitted during functional 
minimization, but r(S) and h(S) are still free. This im- 
plies the additional boundary conditions 0, 0] 



The shape equations are the associated Hamilton equa- 



dr 
dh 



p r (S) and 



s=S 



= Ph{S) 



(12a) 
(12b) 
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The Hamilton equation ||9fJ shows that ph is an integral 
of "motion", and the boundary condition l|12b|) fixes its 
value to zero. Hence, ph drops out of the problem ev- 
erywhere. The condition on p r can be made very useful 
by a little more thought. Since for a > the angle ip 
converges to zero in an essentially exponential way (see 
Sec.E3, the expression for the Hamiltonian will converge 
toward p r . Thus, the requirement of a flat profile implies 
H = H{S) — > p r {S) = 0, or in other words, the pro- 
file is flat if the Hamiltonian vanishes. Using Eqns. I|7a(l 
and (jHJ , this can be turned into a condition for p r at the 
contact boundary: 



ap r (0)= z Z) {l + 2~az-{a^) 2 ) . (13) 

The only remaining variable for which the contact 
value is not yet known is p^, or alternatively ip. At this 
point a bit of care is required. It is well known that the 
balance between adhesion energy and elastic membrane 
deformation results in a boundary condition on the con- 
tact curvature [l^- For curved substrates this becomes 

m 

aip :=a ip(0) = 1 - yfOi . (14) 

However, it is crucial to understand that this condition 
only holds for the final equilibrium shape of the complex. 
In the present case the situation is different, because the 
aim is to calculate the energy E[ lee at any given value 
of the penetration z. Therefore the adhesion balance is 
restricted and Eqn. i|14|) generally will not hold. Never- 
theless, by later imposing dE/dz = for identifying the 
equilibrium degree of penetration, Eqn. (jl4|l is recovered. 
In fact, this would be one way to derive it. 

Rather than Eqn. (|14f> . it is the condition of asymp- 
totic flatness that will determine ipQ. hi practice, this 
can be done via a shooting method: For a trial value of 
ipo integrate the profile, find the arc-length so at which 
the angle vanishes, i. e. ip(so) — 0, and observe at which 
radial distance v(sq) this happens. Now adjust ipo itera- 
tively and search for the value(s) at which r(so) diverges. 
One thereby finds the contact curvature as a function of 
penetration, ipa(z), in what amounts to a nonlinear eigen- 
value problem. 

At this point a side-note seems appropriate. It turns 
out that at a = <r c = 4.721139 . . . the nature of this eigen- 
value problem changes qualitatively, since a region of z 
values emerges (at z c = 1.86289 . . .) for which there ex- 
ist three contact curvatures which yield profiles satisfying 
all boundary conditions (the function ipo(z) develops an 
S-shape) . The correct solution has to be identified based 
on the criterion of lowest energy. For somewhat larger 
values of a the lowest energy solution cannot even be re- 
alized physically: The corresponding curvature ipo m ay 
become larger than 1/a, which is geometrically impossi- 
ble because the membrane cannot bend into the colloid 
it adheres to. This, however, is not a problem in the 
present case of an adhesion balance, since the curvature 
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FIG. 2: Structural wrapping phase diagram in the plane of 
reduced adhesion constant w and reduced lateral tension a, 
in the low tension regime a < 1 close to the triple point "T" 
(w = 4, a = 0). The dashed line "W" marks the continuous 
transition at which partial wrapping sets in, the bold solid line 
"E" indicates the discontinuous transition between partially 
wrapped and fully enveloped, and the short dashed lines "Si" 
and "S2" are the spinodals belonging to "E". The fine dotted 
line w = 4 + 25" close to "E" indicates where the fully wrapped 
state has zero energy. 

boundary condition from Eqn. (|14|l implies that for any 
wrapping geometry in which the point of detachment is 
not fixed the contact curvature has to be smaller than 
1/a. In fact, it turns out that the accessible range of 
multivalued contact curvatures always lies inbetween the 
transition from partially to fully wrapped, hence it has 
no direct consequences on the phase diagram [35]. But 
the mathematical properties of this bifurcation may per- 
mit some insight into the general nature of the solution, 
which, however, will not be pursued in the present paper. 



C. Structural phase diagram 

Numerically performing the calculations indicated 
above yields the shape profile and hence, via Eqn. JSJ, 
the free membrane energy for any given value of a and 
z. From Eqn. Q one then determines the total energy 
as a function of a, w, and z. The minimum in z within 
the range [0; 2] corresponds to the equilibrium state and 
one obtains the phase diagram as depicted in Fig. [21 25] . 

The transition from free to partially wrapped (see 
dashed line "W") remains unchanged compared to the 
case where -Ef ree was neglected. The reason for this is 
that - just like the tension term - the energy of the free 
membrane is of higher than linear order for small z (this 
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FIG. 3: Height of the energy barrier -Ebarricr = £barricr/7rft as 
a function of reduced tension in the low tension regime a < 1, 
at the value of w where the transition from partial wrapping 
to full envelopment occurs. The inset illustrates the shape of 
the function E(z) and defines the concomitant energy barrier 
(<r = 0.5 in this example). 

again follows analytically from a small gradient expan- 
sion, see Sec. IIV|I . The physics is thus determined by a 
balance between bending and adhesion alone. However, 
the transition from partial wrapping to full envelopment 
(the solid line "E" ) changes significantly: As can be seen 
in Fig. 13 an energy barrier separates the fully and par- 
tially wrapped states, rendering the transition discontin- 
uous. This energy barrier turns out to be mostly ten- 
sion (not bending) energy stored in the free membrane 
of partially wrapped colloids, and its height can be quite 
substantial. For instance, with a = 0.02dyn/cm (a typ- 
ical value for a cellular tension 26] ) , a = 30 nm (capsid 
radius of Semliki Forest Virus, as an example for a col- 
loidal particle to be wrapped) and re = 20 fceT, one finds 
a w 0.22 and from that ^barrier ~ 22k B T. This bar- 
rier is too large to be overcome by thermal fluctuations 
alone 36] . However, upon increasing the adhesion en- 
ergy w and thereby going deeper into the region of full 
envelopment, the energy barrier separating the partially 
and the fully wrapped state decreases, ultimately vanish- 
ing at the spinodal line "Si". Conversely, once the col- 
loid is fully wrapped, the same energy barrier prevents 
the unwrapping transition, and one has to decrease the 
value of w further in order to remove this barrier — see 
the second spinodal "S2" in Fig. |2 Cycling across the 
envelopment transition "E" thus gives rise to hysteresis, 
as is illustrated for the particular case a = 1 in Fig. 
for which the energy barrier is -Ebarricr ~ 66 k^T using 
the same system properties as above. Interestingly, this 
hysteresis is so pronounced that one "skips" entirely the 
partially wrapped region upon unbinding. 

Both w and a are proportional to a 2 ; therefore, a scan 
of the particle radius a at fixed values of re, a, and w 
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FIG. 4: Hysteresis loop of crossing the envelopment transi- 
tion by changing the adhesion energy w for the particular 
tension a = 1. At w = 4 binding sets in, at w « 6.1 the fully 
enveloped state becomes stable. However, only at w fa 7.5 
does the huge energy barrier of 66 k^T separating it from 
the partially wrapped state vanish. On the unbinding branch 
the system again remains metastable beyond the actual tran- 
sition, the stable partially wrapped branch of low bindings 
is entirely skipped and "replaced" by a metastable partially 
wrapped branch below w — 4, featuring large values of z. At 
W « 2.7 the energy barrier for unbinding vanishes. 



yields lines in the phase diagram which pass through 
the origin. The shape of the envelopment boundary "E" 
then implies that for w/a < 1.37 particles will not be- 
come fully enveloped, irrespective of their size, while for 
w/a > 2 all sufficiently large particles are enveloped. In 
the small region inbetween, 1.37 < w/a < 2, particles 
are only enveloped if they are neither too small nor too 
large. The asymptotic envelopment condition for small a 
coincides with the boundary at which wrap ping sets in, 
which is a = ^2n/w or w/a = 2 (A/a) 2 [HHg. A * the 
onset of the possibility of full envelopment, w/a ~ 1.37, 
the first particles to be enveloped have a radius a ~ 4.4A. 
Fig- El summarizes these results. 

The energy of the free membrane, Ef tee , vanishes not 
only in the limit z — > but also in the full- wrapping limit 
z — > 2. This is reminiscent of the case of an ideal neck 
connecting two vesicles [33. The reason is essentially 
that the membrane shape locally approaches a catenoid 
[38| . This fact is very convenient, because knowledge of 
the exact energy of the enveloped state greatly simpli- 
fies the discussion of the structural transitions. This will 
be further exploited in Sec. IIVI As an immediate con- 
sequence, it becomes possible to estimate the point of 
envelopment by comparing the known energy of the fully 
enveloped state not with the partially wrapped state but 
simply with the free state, namely, E = 0. This gives 
the boundary w = 4 + 2a, which is also plotted in Fig. [21 
and which actually becomes asymptotic to the real phase 
boundary in the limit a — » (in a complicated logarith- 
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FIG. 5: Influence of the particle radius a on the wrapping 
behavior. Sufficiently large particles will always at least par- 
tially wrap. In the range 1.37 < w/a < 2 particles only 
become enveloped if they are neither to small nor too large. 

mic fashion). Note that this line differs from the phase 
boundary of the case where -Ef re e had been neglected by 
a factor 2 in the slope (i e., the prefactor of a) — and in a 
maybe unexpected way: The region in the phase diagram 
belonging to fully enveloped states grows at the expense 
of partially wrapped states. Even though bending and 
tension energy work against adhesion, they can actually 
•promote wrapping. The reason is that partially wrapped 
states with a large penetration can lower their energy by 
completing the wrapping, which provides another means 
to understand why the transition is discontinuous |39| . 
The same has been found for colloids adhering to quasi- 
spherical vesicles [l6|. 

For increasing a the bending energy should ultimately 
become negligible compared to the tension. Indeed, in 
the limit k — > the term i?f rco vanishes, because the 
membrane is flat immediately after detaching (smooth- 
ness of the slope is no longer required). The equilibrium 
penetration of partially wrapped colloids, as deduced 
from Eqn. is then z = w/2a. This equation can 
be rewritten as w = a[l + cos(7r — a)] and is thereby rec- 
ognized as the Young-Dupre equation [40| , which relates 
adhesion and tension to the contact angle, here ir — a. 
The envelopment boundary "E" consequently occurs at 
w = 4(7, i. e., where the penetration is z = 2, or, equiv- 
alently, where the contact angle vanishes and the mem- 
brane completely "wets" the colloid. 

On the basis of these results one expects a cross-over 
from the small gradient asymptotic phase boundary w = 
4 + 2a, valid below a ~ 1, to a large tension limit w — 
4cr ~ 4 + 4<j. Fig. [fj] confirms this. However, it is quite 
remarkable how many orders of magnitude of variation 
of the reduced tension it takes to establish the transition 
toward the high tension asymptotic: At a — 1 the curve 
is about 0.142 away from the zero tension asymptotic; 
getting as close as that to the high tension asymptotics 




FIG. 6: Phase boundary between partially and fully wrapped 
state plotted against ten decades of reduced tension, a. The 
combination (w — 4)/er is seen to cross over from the value 2 at 
small a, which follows from the small gradient expansion, to 
the Young-Dupre-limit 4, which results when the energy Ef TBe 
of the free membrane is negligible. The two dotted curves at 
small and large a correspond to the small gradient estimate 
(12711 and the scaling prediction 1351 1 . respectively. 



requires a ~ 2 x 10 5 . Over the intervening five orders of 
magnitude of reduced tension, the influence of bending 
and tension cannot be easily disentangled. In Sec. it 
will be shown that the Young-Dupre-limit is reached in 
a power law fashion with a small exponent 1/3, which 
partly explains this slow crossover. 

It may be worth pointing out that the large range of 
values of a is not experimentally unreasonable, because 
each of the three variables entering the reduced tension 
can vary by a few orders of magnitude: Membrane ten- 
sions between 0.01 dyn/cm and 1 dyn/cm are typical 26], 
as are bending constants between 1 k-g,T and 100 k^T |4lj . 
Assuming colloidal radii between 20 nm and 2 /iin yields 
a range for a from about 10 -2 up to 10 6 . 



IV. SMALL GRADIENT EXPANSION 

One particular result from the numerical solution of 
the nonlinear problem is the following: For sufficiently 
small tension the equilibrium penetration shortly before 
envelopment ensues is quite small (see the inset in Fig-EJ) , 
as is the concomitant perturbation of the flat membrane. 
Therefore, this region of the phase diagram should be 
amenable to an approximate treatment of the differential 
equations which corresponds to a lowest order expansion 
around the flat profile. 
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A. Functional and linear shape equations 



If the shape of the membrane is only weakly perturbed, 
a Monge representation giving the profile height ft as a 
function of the position r = (x, y) in the reference plane 
is applicable. Bending plus tension energy can then be 
written as 



E = / d 2 r y/1 + (Vh) 2 



V • 



(15) 

where V is the two-dimensional nabla operator in the 
reference plane. Expanding the two terms in the inte- 
grand up to lowest order in Vft gives the small gradient 
expansion of the energy functional 



E = |d 2 r{|(V 2 /,) 2 + |(Vft) 2 } 



(16) 



The functional variation SE = finally yields the linear 
shape equation [42|: 

V 2 (V 2 - A~ 2 ) h = , (17) 

where A is the length introduced in Eqn. P)l. 

B. Equilibrium profile and energy 




FIG. 7: Exact membrane profiles (solid curves) and small 
gradient approximation (dashed curves) for two fixed detach- 
ment angles a = 30° and a — 60°. The reduced tension is 
cr = 0.1. 



profile is important if one attempts to infer physical prop- 
erties of the membrane or the complex by measuring the 
membrane deformation and working backwards. Using 
the linearized prediction of the profile then would lead to 
incorrect conclusions, for instance to an underestimation 
of the degree of wrapping. 

The energy corresponding to the optimized membrane 
shape is obtained by inserting the profile (ftT?)) back into 
the functional (|16fl . The necessary integration can be 
performed analytically, yielding [25| 



Ehc 



k 3 \ K a (ka/X) 



1 - k 2 J Ki (ka/X) 



(20) 



The differential equation H17|) is solved by eigenfunc- 
tions of the Laplacian corresponding to the eigenvalues 
and A -2 . In the present cylindrical symmetry the general 
solution can therefore be written as 

h(r) = h 1 +h 2 ln(r/A) + h 3 K {r/X) + h 4 I Q (r/\) , (18) 

where Kq and Io are modified Bessel functions. Since 
loir) diverges as r — ► oo, the condition of a flat profile 
requires hi = 0. And the coefficient hi has to vanish 
since otherwise the energy density is not integrable for 
r — > oo. The two remaining coefficients are obtained by 
fixing height and slope of the profile at the point where 
it touches the colloid. A straightforward calculation then 
gives the small gradient profile 

Kr) 1 | A k K Q (ka/X)-K Q (r/X) 

a a 1 — z Ki(ka/X) 

where the abbreviation k = sin a = y z(2 — z) has been 
used. Fig. \7\ gives an example how the small gradient 
prediction of the profile compares to the full solution. If 
the detachment angle a is sufficiently small, the overall 
membrane deformation remains also small, and the pro- 
file from the linearized theory follows the full solution 
quite accurately. However, for a somewhat larger a sig- 
nificant deviations appear: the membrane deformation is 
predicted to be substantially larger than it actually is. It 
is worth pointing out that a good understanding of the 



C. Exact asymptotic results 

The small gradient expansion becomes asymptotically 
exact in the limit of vanishing penetration, z — » 0. It is 
then also permissible to replace the expression H20|l by 
its small z expansion 

E ircc = -2~oz 2 (2 7 + ^ y) + °( z3 ) ' ( 21 ) 

where 7 = 0.5772 ... is the Euler-Mascheroni constant. 
Note that up to an (important) logarithmic correction 
this has the form of a tension energy. Inserting this into 
Eqn. (0} one obtains the small gradient free energy, up to 
quadratic order in the penetration z, and one can discuss 
the phase behavior. 

As already mentioned in Sec. IIII CI the energy of the 
free membrane is of higher than linear order for small 
z, hence the transition from free to partially wrapped is 
dictated by the balance between bending and adhesion 
alone, giving the phase boundary w — 4. The equilib- 
rium penetration follows from dE/dz — 0, which can be 
rewritten as 

We w^_^± e 2 7 with M/ . = 27 + ln ^£^ (22 ) 
8 2 

The solution of this transcendental equation is known as 
the Lambert W function |43( , and in the present situation 
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it is the branch —1 which is needed. One thus obtains 
~- 2e " 27 ^- 2e_2 Ve- = -^i (23) 



with 



aW 



W = W(w) = W- 



4aW 



w — 4 



,20 



(24) 



Equation 11'- ill thus gives the penetration as a function of 
w and a. For x — > 0~ the function W_i(a;) diverges to 
— oo in a logarithmic way |44|. therefore the penetration 
z(w) increases at w = 4 for all values of the tension 
linearly up to a logarithmic correction. Note also that 
its dependence on the tension is very simple, namely, 
just inversely proportional. 

The expression (|2*3|) is of course not valid for all w > 4, 
because the system must also cross the transition toward 
the enveloped state. At this point one has to make use of 
a piece of information known from the nonlinear studies, 
namely, the energy of the fully enveloped state. It was 
found that the energy of the free membrane vanishes as 
z approaches 2, hence the full envelopment boundary is 
given by the simultaneous solution of dE/dz = and the 
additional equation 



E(z) = E{2) = -2(w - 4) + 4a 



(25) 



By eliminating the logarithmic term between those two 
equations one obtains a quadratic equation for z. After 
inserting its solution into Eqn. (|23[1 . the final expres- 
sion can be solved for a, and one arrives at the phase 
boundary 



w - 4 




1 



2W {2Wf 



(26) 



Remembering the divergence of W as w — > 4 + , Eqn. H26|) 
shows that in the limit of weak binding the phase bound- 
ary approaches w — 4 + 2a, as has been anticipated from 
the numerical results in Sec. IIII CI (see also Fig. |2J. Us- 
ing the lowest order approximation W_i(x) ~ hi\x\ at 
x — > 0~ and expanding the square root in Eqn. (|26|) . 
on gets the approximate asymptotic phase boundary 



w - 4 



1/ , w - 4\~ 1 
l + -(27 + ln— - 



(27) 



which is indicated by the left dotted curve in Fig.[fJ] The 
expression l|26|l is significantly more accurate, but it re- 
quires the function W-i to be evaluated. In any case, 
since W_i(x) is only real for — 1/e < x < 0, even the full 
expression (|2*6*)l exists only up to w — 4+8 e -1-27 « 4.928 
or, equivalently, a = e" 1 " 27 /(2 - y/3) » 0.433. At the 
upper boundary for a one finds z = 2(2 - V3) « 0.536. 
Larger penetrations than this cannot be described within 
the quadratic approximation l|21|) to the small gradient 
energy 

Finally, the value of the penetration z on the discon- 
tinuous phase boundary can be obtained by eliminating 



w between the two defining equations dE/dz — and 
(12511 . Solving the remaining equation for a yields 



a = — exp 

z 



4 - z 1 n \ z«i 2 



(28) 

The second approximate relation can also be solved in 
terms of the Lambert W function: 



CT<1 



1 



2W_i(-er/4) 





-In 






H 















(29) 

In the limit a — > 0, i. e., when approaching the triple 
point, the penetration on the discontinuous phase bound- 
ary vanishes. Hence, the jump in order parameter ap- 
proaches 2, i. e., the transition becomes increasingly dis- 
continuous at smaller a. However, Fig. |3 demonstrates 
that the barrier vanishes in the limit a — > 0, so from 
this point of view the transition becomes more continu- 
ous. The triple point (w = 4; a — 0) is thus quite un- 
usual. Another peculiarity is that along the phase bound- 
ary z does not approach in an algebraic way; rather, 
Eqn. (|28|l shows that z(a) has an essential singularity 
at a = 0. All this is related to the fact that the small 
gradient expression for the energy is not a conventional 
Landau expansion in powers of the order parameter, z, 
since the quadratic term has an additional logarithmic 
factor. This lies at the heart of all logarithmic correc- 
tions encountered above (manifest also in the occurrence 
of the function W), and it renders the standard classifi- 
cation schemes for critical points inapplicable here. 

The exact asymptotic phase boundary can be obtained, 
because information about the energy of the fully wrap- 
ped state is available. However, for the barrier the sit- 
uation is different: Even if the equilibrium penetration 
is very small, the location of the barrier (i e., the pene- 
tration Zbarrior at which the energy is largest) occurs at 
large z (see e.g. the inset of Fig. 0). In fact, numerical 
evidence suggests that lim^o ^barrier = 1 from above. It 
is therefore impossible to obtain the height of the barrier 
by extending the above small gradient analysis. 



V. SCALING LAWS IN THE HIGH TENSION 
LIMIT 

As the tension grows, so does the equilibrium penetra- 
tion on the phase boundary toward the fully enveloped 
state. The location of this transition can then no longer 
be obtained within the small gradient framework of the 
previous section. Still, the numerical results in Sec. IIII CI 
strongly suggested that the system displays a well defined 
and simple asymptotic behavior in the high tension limit 
(see, e.g., Fig. EJ. Unfortunately, treating the curvature 
as a small perturbation to the tension is tricky, because 
this leads to a so-called "boundary layer problem" : The 
solution features a finite variation over a range which van- 
ishes in the perturbative limit |45| . In the present case, 
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the membrane has to bend away from the colloid toward 
the flat plane (i. e., ip(s) has to change from a to 0) in 
a region of vanishing arc-length. Typically, such prob- 
lems are dealt with by a subtle matching procedure (an 
example is provided by the treatment of an ideal neck in 
Ref. |23). Somewhat less ambitious, the current section 
shows how the asymptotic behavior can be quantified by 
starting with reasonable scaling assumptions about the 
boundary layer. Still, the resulting formulas will turn out 
to be remarkably robust. 

A useful observation to start with is that for large a the 
equilibrium penetration approaches z = w/2a (see the 
discussion of the Young-Dupre-limit k — ► in Sec. lIII 
Using the contact curvature boundary condition Q14JI. 
this would predict the asymptotic relation 



aipo = 1 — \AB 



- <r»i 
w ~ 



(30) 



This suspicion is indeed confirmed by a check with the 
numerical results (data not shown). The proportional- 
ity to \fd = a/A is not surprising, since A is the typical 
length on which the membrane bends. However, the pro- 
portionality to \[2z = 2 sin ^ is not obvious jig . 

Equation l|30|) can be used to infer the asymptotic be- 
havior of several more variables, by virtue of the follow- 
ing scaling argument. Its aim is to estimate the energy 
of the free part of the membrane, which for large a is 
largely stored in a small toroidal rim at contact (this is 
the boundary layer). This toroid has the axial radius 
a sin a and a typical meridinal radius which scales like 
1/t/V Its area is thus proportional to (asxn.a)/ipo, and 
the tension contribution becomes 



*C ~ o x 



(31) 



The two principal curvatures are 1/et and ipo , where the 
second one clearly dominates in the high tension limit. 
Hence, the bending energy of this torus scales like 



TP 



bond 
free 



a sin a 



KX - 



x (ipo) 2 ~ kV^z^/2 — z ~ K\fd\J 2 — z 

where in the last step the prefactor z has been dropped, 
since for high tension the equilibrium penetration at the 
transition is close to 2. Eqns. (|31(l and (|32|l show that 
in the limit of large tension and close to full wrapping 
the energy of the free membrane can be written in the 
following way 



E flee ~ 2A\f& V2~ 



(o>M«2) 



(33) 



where the proportionality factor A does not depend on 
a or z, and the additional factor 2 is included for conve- 
nience. 

It is worth pointing out that the z-dependence of the 
scaling form in Eqn. I|33|l can also be understood in the 
following intuitive way: If one conceives of the strongly 
curved region at detachment as giving rise to a line en- 
ergy, it follows that Ef vcc ought to be proportional to 
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FIG. 8: Scaling plot for the penetration z on the phase bound- 
ary (solid) and the location of the barrier, Zbarricr (dashed). 
According to Eqns. I).'-{4|I and i'M>l the combination (2 — z)a x ^ 
should approach a constant value, which should be the same 
for both cases if the latter is multiplied by the additional fac- 
tor 2/(2 — y/S). The dotted line indicates the asymptotic limit 
A 2/3 « 3.17. 



the length of this line, which is a-J 'z(2 — z). For z close 
to 2 this has the same characteristic variation y2 — z 
as Eqn. I|33|) . However, it must be noted that this form 
holds only in the double limit of large tension and large 
penetration. Generally, -Ef rcc is not well represented by 
a simple line energy alone. 

One can now insert the expression l)33|l into Eqn. Q 
and discuss the phase behavior. Eliminating w between 
the two equations E(z) = E(2) and dE(z)/dz — gives 
the penetration z at the transition as a function of re- 
duced tension: 



Z** 1 2 - J 42/3~-l/3 



(34) 



The high tension limit of the penetration is thus reached 
in an algebraic way with an exponent —1/3; see Fig. [S] 
Eliminating z instead of w gives the envelopment bound- 
ary: 



w — 4 5>1 



4-3A 2/3 ct" 1/3 



(35) 



showing that its asymptotic value is also reached alge- 
braically with an exponent —1/3; see Fig. El 

Information on the barrier can be obtained by further 
studying the scaling energy. After inserting the phase 
boundary (|35|l back into the energy one determines the 
location of the maximum via dE(z)/dz = 0: 



Zba 



V3 



A 2/3~-l/3 



(36) 



The location of the barrier thus reaches the asymptotic 
value 2 in the same way as the location of the transition, 
only the prefactor is different by (2 — y/3)/2. Fig.[S]also 
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FIG. 9: Barrier height for the transition between partially 
wrapped and fully enveloped as a function of reduced tension 
a on a double logarithmic scale. The two dotted lines indicate 
the two different regimes: For high tension the barrier height 
scales with an exponent 1/3, while an empirical power-law fit 
to the low tension regime gives the value 0.86 |47l| . 

shows a scaling plot of the location of the barrier, in 
which this additional factor has been explicitly included. 
The fact that both curves in Fig. [S] approach the same 
limit indicates that the present scaling argument predicts 
more than the exponent: it correctly predicts the ratio 
of the prefactors as well. 

Finally, the barrier height is determined as the differ- 
ence between -E(-Zbarricr) and E(2), for which one finds 

B baIrie /i 1 ^(2V3-3)A 4 /3^/3. (37) 

This is illustrated in Fig. OH The two asymptotic power 
laws meet at the crossover point CT cr oss ~ 4.72. The nu- 
merical value is intriguingly close to the critical tension 
<7 C mentioned at the end of Sec. lIIIBl but this is probably 
coincidental. 

That the above scaling argument gives the correct re- 
lation between the prefactors can also be checked in the 
following way: Each of the equations 1)34(1 . i|35|) . (|36[1 . and 
(|37|l describes a scaling relation for a different variable, 
but the prefactors all involve A. Numerically one can 
determine A by an asymptotic fit to the high tension val- 
ues of these four variables, determined from the nonlinear 
studies of Sec. IIIII In all cases one finds the same result: 
A w 5.650. The scatter among the four results relative 
to the average value is very small, only about 6 x 1CP 4 . 



VI. DISCUSSION OF A BIOLOGICAL 
EXAMPLE 

In the previous three Sections a theoretical description 
of the adhesion and wrapping behavior between a colloid 



and a fluid membrane in terms of continuum elasticity 
theory has been developed, using the full nonlinear shape 
equations, their small gradient expansion, and scaling ar- 
guments. In this final section the results obtained are 
used to again make contact with a biological application 
of such wrapping events mentioned in the Introduction, 
namely, the maturation of animal viruses by budding. 

While the reduced tension a used throughout the paper 
can span many orders of magnitude, it is important to 
realize that in a biological context the variation is more 
restricted. Tensions of cellular membranes reported in 
the literature vary between 0.003 dyn/cm and ldyn/cm 
|26j . Larger values soon result in a structural failure of 
the bilayer. On the other hand, typical bending con- 
stants of membranes are in the range of a few tens up 
to about a hundred k^T. From these numbers one finds 
that the characteristic membrane length A from Eqn. 
varies roughly between ten and at most a few hundred 
nanometers. Interestingly, this coincides with the range 
of particle sizes for which the scenario treated in this pa- 
per is biologically meaningful, for quite different reasons: 
Particles much smaller than 10 nm are more likely to be 
transported across a biomcmbrane by means of chan- 
nels, while wrapping of particles much bigger than a few 
hundred nanometers can no longer be described without 
considering the concomitant significant rearrangements 
of the cytoskeleton. Hence, if wrapping events of the 
kind discussed in this paper take place on cellular mem- 
branes, they are bound to occur in the regime in which 
the reduced tension a = (a/A) 2 is of order 1. 

A prominent class of colloidal particles exactly within 
the right range, for which such wrapping events occur and 
have been studied in great detail, are the nucleoprotein 
capsids of many animal viruses belonging for instance 
to the families of Togaviridae, Coronaviridae, Retroviri- 
dae, Rhabdoviridae, Ortho- and Paramyxoviridae, and 
Hepadnaviridae 2]. During their final maturation step 
the viral capsids are enveloped by a cellular membrane 
{e.g. the plasma membrane or the endoplasmic reticu- 
lum) in an event which is believed to be independent of 
active cell processes and by which they ultimately leave 
their host. In the simplest case adhesion is due to a direct 
interaction between the capsid and the membrane (for in- 
stance in the case of type D retroviruses [2]). However, 
more common is an adhesion mediated by viral trans- 
membrane proteins (usually called "spikes") which can 
attach at specific binding sites on the capsid 0, 0, Eli 
and for which Scmliki Forest Virus (SFV) is the classical 
example. SFV has a capsid radius of about 30 nm and 80 
spikes. Assuming a typical membrane bending stiffness 
of k ~ 20 k^T, one finds that the wrapping boundary 
"W" at w = 4 corresponds to a binding (free) energy per 
spike of about 6 /cb? 1 , which is physically reasonable. 

The above estimates indicate that viral wrapping 
events can be expected to take place in the low ten- 
sion regime of the phase diagram, close to the two phase 
boundaries. This pattern is found time and again in bi- 
ology: Systems often seem to have evolved to lie close to 
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phase boundaries, because this permits large "effects" to 
be triggered by comparatively small parameter changes. 
Recall, however, that the envelopment transition has 
been found to be discontinuous and associated with a 
substantial energy barrier, which nature somehow has to 
overcome. A conceivable solution of this problem would 
be provided by a coupling between curvature and com- 
positional degrees of freedom |48|, |49|, |50|, |5JJ with the 
result of enhancing the concentration of lipid species in 
the highly curved rim which actually prefers a high cur- 
vature. This would lower the rim energy and thereby the 
wrapping barrier. 

It is crucial to realize that it is biologically feasible 
to actually move in the phase diagram of Fig. For 
instance, cells actively control and adjust their surface 
tension for the purpose of surface area regulation [26(. 
Even more dramatic changes in tension can occur when 
one switches between adhering membranes. If viral cap- 
sids get spontaneously wrapped, they evidently must be 
in a region of the phase diagram in which the wrapped 
state is stable (and, moreover, in which it is not rendered 
inaccessible by a large barrier). But the virus cannot stay 
wrapped forever. As it infects a new host cell, it typi- 
cally becomes internalized via receptor mediated endocy- 
tosis and ends up in an endosome, which it again has to 
leave in order to avoid ultimately being digested by cel- 
lular toxins. Many viruses leave the endosome by fusing 
their outer envelope with the endosomal bilayer, but if 
the capsid were too strongly attached to the membrane, 
it could not be freed this way. It is usually assumed 
that the biochemical changes within the endosome which 
lead to the fusion event (in particular, a lowering of pH) 
also diminish the strength of adhesion. However, within 
the theoretical framework established in this paper it is 
tempting to speculate about an alternative mechanism: 
If the bilayer tension of the endosome is larger than the 
tension of the membrane at which the capsid became en- 
veloped, unwrapping can be efficiently promoted by ver- 
tically crossing the phase boundary "E" from enveloped 
to partially wrapped, as can be seen in Fig. More- 
over, the horizontal adhesion axis of the phase diagram 
can not only be changed by chemically modifying the 



spikes, but also by controlling their density in the mem- 
brane |52| . This may not only be relevant in the initial 
wrapping event, in which an increasing density of spikes 
in the membrane can push the system over the envel- 
opment boundary, but also in the unwrapping process, 
when after fusion the spikes can readily diffuse into the 
essentially spike-free endosomal membrane and thus re- 
duce the binding free energy. 

The above example illustrates how the physical prin- 
ciples discussed in this paper can be directly relevant 
in a biological context. Unfortunately it is often hard 
to disentangle them from other biological processes or 
secondary effects of the experimental set-up. Hence, a 
quantitative test of the present work appears more prac- 
ticable in well controlled model systems, e. g., similar to 
the ones studied in Refs. 0,0]. Nevertheless, the phys- 
ical results presented here can provide valuable insight 
into biological problems which may complement other 
approaches. As an example one might think of a way to 
measuring cellular tension which is an alternative to the 
current method of pulling a tether 0, [2^. The above 
analysis has shown how the degree of wrapping of a col- 
loid depends on the applied tension — in the regime ac- 
cessible by the small gradient expansion it is simply in- 
versely proportional, see Eqn. J22Jl. One can thus use 
suitably coated colloids as tension probes. Unlike the 
tether approach this method is in principle also applica- 
ble to intracellular membranes, even though a noninva- 
sive determination of the degree of wrapping will be very 
difficult for small beads. The theory developed in this 
work should then be useful for analyzing the results of 
such measurements. 
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